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The electronic and optical properties of (Zn,Mg)(S,Se) wide-gap solid solutions are studied using 
ab initio techniques and starting from the previously determined atomistic structure of the alloy. 
Compositional disorder is shown to close substantially the gap with respect to the predictions of the 
virtual-crystal approximation. The bowing of the fundamental gap vs. composition predicted by our 
calculations is in very good agreement with experiments available for the Zn(S,Se) pseudo-binary 
alloy. At temperatures typical for MBE growth, the quaternary alloy displays a rather large amount 
of short-range order whose effect is to slightly but unmistakably open the gap. Our results agree 
well with recent experimental data for the quaternary alloy. 



In the last decade much effort has been devoted to the 
study of semiconducting materials capable of operating 
in the short wavelength optical range. The main goal 
of these efforts is the realization of light-emitting and 
laser diodes encompassing the entire visible-light win- 
dow and, in perspective, the industrial-scale production 
of high-density storage optical disks and light sources for 
full-color displays. II- VI and nitride semiconductors are 
the most promising among these materials. Among II- 
VI semiconductors, ZnSe-based materials have a special 
importance (IJ], and quaternary Zn^Mg-^^SySei-j, alloys 
are commonly used as cladding layers in II- VI blue-green 
laser diodes. This system — which was introduced a few 
years ago — has the desirable property that its lattice 
parameter, ao, and its fundamental band gap, E g , can 
be tuned fairly independently by acting on the concen- 
trations x and y. In spite of the technological impor- 
tance of this material, extensive experimental studies of 
its electronic and optical properties over a wide range of 
compositions and of their dependence on the structural 
and thermodynamic equilibrium properties are still lack- 
ing. 

We recently reported (3) on a theoretical study of the 
thermodynamical properties of this system, performed 
using state-of-the-art ab-initio methods. In that study 
the phase diagram of the quaternary alloy was deter- 
mined, and the homogeneous alloy was found to be stable 
against segregation or the formation of ordered struc- 
tures at temperatures typical of molccular-beam epitaxy 
(MBE) growth (« 550 K). Nevertheless, a large amount 
of short-range order (SRO) characterized by the occur- 
rence of Zn-Se and Mg-S clustering among first-nearest 
neighbors was found to occur and to persist even at very 
high temperatures (~ 1700 K). 

In this Letter we report on the first extensive theo- 
retical study of the electronic and optical properties of 
(Zn,Mg)(S,Se) alloys both in the case of the experimen- 



tally well studied pseudo-binary ZnS^Sei-a; system, and 
for the quaternary solid solution in conditions of lattice- 
matching to GaAs substrates. In the latter case the pres- 
ence of the short range correlations will be shown to in- 
duce a small but significant opening of the fundamental 
band gap. 

Our work is based on density-functional theory (DFT) 
within the local density approximation (LDA), and the 
plane-wave pseudopotential method. All the technical 
details of our calculations are the same as in Ref. [[|. 
DFT-LDA is known to underestimate electronic excita- 
tion energies, and many-body corrections should be con- 
sidered, e.g. using the GW method Q. In spite of this, 
the agreement between the shape of the band-structure 
predicted by the LDA and that measured experimentally 
is usually satisfactory, and semi-quantitative agreement 
can be simply obtained by a rigid upward shift of the 
conduction bands, using the so-called scissor operator. 

Ideally, the study of the electronic properties of an 
alloy should proceed in three steps: i) the determina- 
tion of the equilibrium structural properties of the sys- 
tem, including the relevant atomic correlations, as we 
obtained and discussed in Ref. ||; ii) the calculation of 
the alloy band structure, taking into account the atom- 
istic structure previously determined; Hi) the inclusion 
of electronic correlations. It is still very hard to combine 
many-body calculations Q with a proper description of 
disorder, and in the following we assume that disorder 
and quasi-particle effects can be treated independently: 
we deal explicitly with compositional disorder at the LDA 
level, and we add semi-empirically the many-body correc- 
tions by a simple interpolation of the scissor operators 
appropriate for the pure compounds. 

Starting from the structural data determined in Ref. 
, we have calculated the electronic properties of the 
alloys using the special quasi-random structures ap- 
proach H . The latter is based on the observation that — 
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FIG. 1. Energy gap of the Zn(S,Se) pseudo-binary alloy 
as a function of sulphur concentration, x. Solid line: ex- 
perimental data from Ref. [7]. Open diamonds: our theoret- 
ical results obtained neglecting many-body effects. Dashed 
line: quadratic interpolation to our theoretical results. Filled 
diamonds: theoretical results semi-empirically corrected for 
quasi-particle effects (see text). 

for any given composition — atomic disorder mainly af- 
fects the electronic properties of an alloy through the 
short-range atomic structure. A disordered solid solution 
can thus be mimicked by using reasonably small super- 
cells that reproduce the alloy SRO in the first few (typ- 
ically, four or five) shells of neighbors. We have general- 
ized the original method so as to account for any given 
short-range order 0, arbitrary compositions and dou- 
ble sub-lattice disorder. We have verified that atomic 
correlations beyond the fourth shell of neighbors affect 
only negligibly the band structure of the alloys presently 
studied. The supercell structures used in our simulations 
were thus obtained by a simulated-annealing procedure 
by which we searched among 64-atom simple-cubic struc- 
tures the ones which give the most similar pair correla- 
tions up to the fourth shell of neighbors, as compared to 
those obtained from Monte Carlo simulations m . 

The electronic properties of the ZnS x Sei- x pseudo- 
binary alloy, which are well characterized experimentally, 
have been studied by the above techniques, at three dif- 
ferent concentrations, i.e. x =0.25, 0.50, 0.75. In Fig. |l| 
we display our results for the fundamental band gap of 
the alloy and compare them with available experimen- 
tal data 0. As previously discussed, LDA calculations 
systematically underestimate the optical gap. Adding 
a simple linear interpolation of the errors done at the 
pure-compound extremes, brings our theoretical predic- 
tions in very good agreement with experiments. Upon 
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in Ref. [3] for T= 550 K 



this semi-empirical correction, our data can be accurately 
described by the formula: 



E g (x)=E^(l 
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where the curvature, commonly known as bowing param- 
eter, is b ~ 0.56 eV. Our result is in good agreement 
with experiments [uT^], according to which 0.41 eV < 
b < 0.68 eV. It is worth noting that calculations done 
in the virtual crystal approximation (i.e. by neglecting 
disorder and averaging the atomic pseudopotentials in 
the anionic sub-lattice) would predict a small and nega- 
tive bowing parameter: b « —0.1 eV. The accuracy by 
which our calculations predict the non-linear dependence 
of the band gap upon composition in the pseudo-binary 
ZnSzSei-^ alloy indicates an almost linear dependence of 
quasi-particle effects on the concentration, which is en- 
couraging for the study of the electronic properties of the 
quaternary alloy, where the available experimental data 
are scarcer. 

The quaternary alloys displays a large amount of short- 



TABLE I. Energy gap vs. compositions in (Zn,Mg)(S,Se). 
'Virtual': theoretical results obtained within the vir- 
tual-crystal approximation. 'Random': super-cell calcula- 
tions for the random alloy. 'SRO': supercell calculations done 
including SRO, but neglecting many-body effects. 'Corrected 
SRO': as above, but corrected for quasi-particle effects (see 
text). 'Expt': experimental data from Ref. [9]. 
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range order. This is illustrated in Fig. 2 where we display 
the atomic correlation functions calculated in Ref. || at a 
temperature of 550 K and for x = y = ~ • The correlation 
function between two atomic species, X and Y, is defined 
as: 



C XY (r) = (C(r)r ? (0))-(C)(r ? ), 



(2) 



where £(r) and rj (r) are integer random variables whose 
value is 1 if the atomic site at position r is occupied by 
an atom of species X or Y respectively, and otherwise; 
(•) indicates the average over disorder; and (£) = x and 
(rj) = y are the concentrations of the X and Y atomic 
species. The large positive value of the Zn-Se and Mg~ 
S first-nearest neighbor correlations indicates a strong 
tendency to form ZnSe and MgS local clusters. This 
kind of clustering acts so as to reduce the elastic en- 
ergy of the system because ZnSe and MgS have lattice 
parameters similar to each other and to that of the al- 
loy (aznSe = 10. 57a. u. and OMgS = 10. 79a. u.), but very 
different from those of ZnS and MgSe (az n s = 10.09a.it. 
and aMgSo = ll.32a.it.). This tendency also shows in 
the second-nearest neighbor correlations which favor like 
cations and anions. 

The electronic structure of the quaternary alloy has 
been studied for three different concentration pairs along 
the line of lattice-matching to GaAs (i.e. for pairs of 
concentrations, x and y, such that the average lattice 
parameter of the alloy is the same as that of GaAs). 

In Table [j] we report the values of the fundamental 
gap as calculated taking into account SRO, neglecting 
SRO, and in the virtual-crystal approximation, and we 
compare these results with recent experimental data || . 
Quasi-particle effects are taken into account by the same 
kind of semi-empirical corrections described in the case 
of the pseudo-binary alloy. In the present quaternary 
case, the values of scissor operators for the four pure 
compounds are bi-linearly interpolated in-between. Our 
calculated values for the energy gaps are in very good 
agreement with experiments (within 0.04 eV). As it was 
the case for the Zn(S,Se) pseudo-binary alloy, the pre- 
dictions of the virtual-crystal approximation are rather 



poor, resulting in too large a gap. The effect of SRO, as 
compared to the perfectly random solution, is to slightly 
open back the gap. The effect is small, but sizeable when 
disorder is maximum (i.e. for x = y = i). The effects 
of SRO on the alloy band structure can be understood 
qualitatively considering the actual alloy as a perturba- 
tion with respect to the appropriate virtual crystal (VC), 
as done by Baldereschi and Maschke |l(| in an early study 
of the Gaa;Ini_ x P alloy. Consider the actual alloy poten- 
tial as a perturbation with respect to the VC potential: 



Vaiicy {r) = V vc (r) + AV(r). 



(3) 



To second order in the perturbation, the alloy band struc- 
ture is given by 

0n , ^ \{n,k\AV\ri ,k')| 2 

where \n, k) are the virtual-crystal eigenvectors with en- 
ergy e„(k), for the appropriate given composition pair 
(x,y). Since occupied (empty) states mostly interact 
with the other occupied (empty) states, that are closer 
in energy, the effect of the perturbation is to push the 
band edges in the gap region, thus reducing the vir- 
tual crystal gap, in agreement with our findings for 
Zn x Mgi_ x S y Sei_y, random and correlated, quaternary 
alloy as well as for ZnS x Sei_ x alloy. In the presence of 
SRO, this result can be generalized by showing that the 
magnitude of the gap reduction depends quadratically on 
the strength of the localized perturbations which trans- 
form the virtual crystal into the actual alloy, multiplied 
by the appropriate atomic correlation functions. From 
Fig. 2 it can be seen that the main difference between 
the correlated and the random alloy is that the former 
has a larger number of bond-length preserving Zn-Se and 
Mg-S nearest-neighbor pairs and a reduced number of 
bond-stretching Zn-S and Mg-Se nearest-neighbor pairs 
relative to the random alloy. The former type of pertur- 
bation of the virtual crystal is probably weaker than the 
latter, thus explaining the smaller closure of the VC gap 
in the correlated alloy as compared to the random case. 

In order to better characterize the optical gap of the 
alloy, we consider the spectral weight: 



A(E, k) = J2\ <^|P(k)|V„) \ 2 S(E - en), 



(5) 



where ip n and e„ are the wave-functions and energy lev- 
els of the supercell, and -P(k) is the projector over the 
states of quasi-momentum k. We have found that the 
alloy wave-functions at the top of the valence and at the 
bottom of the conduction bands actually have a strong 
T character, thus indicating that the direct-gap nature 
of the pure materials is thus conserved also in the al- 
loy. An analysis of the different atomic contributions to 
the DOS has been carried out by projecting the super-cell 



eigenf unctions onto the atomic-like localized orbitals. We 
found that the states at the top of the valence band have 
a strong anionic character (26% centered on S and 67% 
on Se). The anionic character of the top of the valence 
band is a direct consequence of the ionic character of 
the alloy, while the larger contribution of the Se atomic 
orbitals is consistent with Harrison's tight binding ta- 
bles which indicate that the top of the valence band 
states of S-compounds lies about 0.5 eV below the anal- 
ogous states of the corresponding Se-compounds. The 
low-lying conduction states are rather delocalized, with 
contributions from all the four atomic species and a slight 
prevalence of cationic character (28% Zn and 33% Mg). 

A more complete account of this work can be found in 
the Ph.D. thesis of one of us fl2]| . We thank A. Franciosi 
and S. Rubini for prompting our interest in this prob- 
lem and for many fruitful discussions. This work has 
been partly done within the Iniziativa Trasversale Cal- 
colo Parallelo of the INFM. 
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